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EXIT PROBLEM OF A TWO-DIMENSIONAL RISK PROCESS 
FROM THE QUADRANT: EXACT AND ASYMPTOTIC RESULTS 

By Florin Avram, 1,2 Zbigniew Palmowski 2 ' 3 
and martijn r. plstorius 1 ' 4 

Universite de Pau, University of Wroclaw and King's College London 

Consider two insurance companies (or two branches of the same 
company) that divide between them both claims and premia in some 
specified proportions. We model the occurrence of claims according 
to a renewal process. One ruin problem considered is that of the 
corresponding two-dimensional risk process first leaving the positive 
quadrant; another is that of entering the negative quadrant. When 
the claims arrive according to a Poisson process, we obtain a closed 
form expression for the ultimate ruin probability. In the general case, 
we analyze the asymptotics of the ruin probability when the initial 
reserves of both companies tend to infinity under a Cramer light-tail 
assumption on the claim size distribution. 

1. Introduction. 

The multidimensional renewal risk model. In collective risk theory, the 
reserves process X of an insurance company is modeled as 

(1) X(t)=x+pt- S(t), 

where x denotes the initial reserve, p is the premium rate per unit of 
time and S(t) is a stochastic process modeling the amount of cumula- 
tive claims up to time t. Taking S to be a compound Poisson or corn- 
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pound renewal process yields the Cramer-Lundberg model and the Sparre- 
Andersen model, respectively. Recently, several authors have studied ex- 
tensions of classical risk theory toward a multidimensional reserves model 

(1) where X(t),x,p and S(t) are vectors, with possible dependence be- 
tween the components of S(t). Indeed, the assumption of independence 
of risks may easily fail, for example, in the case of reinsurance, when in- 
coming claims have an impact on both insuring companies at the same 
time. In general, one can also consider situations where each claim event 
might induce more than one type of claim in an umbrella policy [see Sundt 
(1999)]. For some recent papers considering dependent risks, see Dhaene 
and Goovaerts (1996, 1997), Goovaerts and Dhaene (1996), Miiller (1997a, 
1997b) and Denuit, Genest and Marceau (1999), Ambagaspitiya (1999), 
Dhaene and Denuit (1999), Hu and Wu (1999) and Chan, Yang and Zhang 
(2003). 

Model and problem. In this paper we consider a particular two-dimensi- 
onal risk model in which two companies split the amount they pay out of 
each claim in proportions 5± and 82 where 8\ + 62 = 1 , and receive premiums 
at rates c\ and C2, respectively. Let £/» denote the risk process of the ith 
company 

Ui{t) := -5iS(t) + at + u u i = l,2, 

where U{ denotes the initial reserve. We will study here the eventual ruin 
probabilities in two cases: 

1. the Levy model, obtained by taking S(t) to be a general Levy process; 

2. the Sparre-Andersen/renewal risk model, where S(t) is 

N(t) 

(2) S(t) = £ a h 

i=l 

N(t) is a renewal process with i.i.d. interarrival times Q, and the claims 
<7j are i.i.d. nonnegative random variables independent of N(t). 

The intersection of the two cases is the classical Cramer-Lundberg model, 
where S(t) is a compound Poisson process with nonnegative jumps. 

We shall denote by F{x) the distribution function of the "claims" <7j, and 
by A and fj, the reciprocals of the means of Q and <7j, respectively. 

As usual in risk theory, we assume that Ui(t) — > 00 a.s. as t — > 00 (z = 1,2). 
In the case of the Sparre- Andersen model, this amounts to pi > p := ^ = 
Ea/EQ. 

We shall assume that the second company, to be called reinsurer, receives 
less premium per amount paid out, that is, 

(3) pi = — > — =p 2 . 

01 02 
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Several ruin problems will be considered here: 

1. The first time r or when (at least) one insurance company is ruined, that 
is, the exit time of (Ui(t), U2(t)) from the positive quadrant 

(4) Tar(ui,U2) := inf{t > : Ui (t) < or U 2 (t)<0}. 

2. The first time r s j m when the insurance companies experience simultaneous 
ruin, that is, the entrance time of (Ui(t), U 2 (t)) into the negative quadrant 

(5) r sim (ni,n 2 ) := inf{t > : U^t) < and U 2 (t) <0}. 

The associated ultimate/perpetual ruin probabilities will be respectively de- 
noted by ^ r(«i,«2) and Vsim(«i,«2) 

lpor(ui,U 2 ) = P(T or (u 1 ,U 2 ) < CO), 
1p sim (ui,U 2 ) = P(T sim (ui,U 2 ) < OO). 

Letting Tj(itj) = inf{t > : Ui(t) < 0}, i = 1, 2, we also will consider 

(6) Vw(/ui,U2) =P{t\{u\) < oo and t 2 {u 2 ) < oo). 

Denoting by tpi(u) := P(rj(n) < oo), the ruin probability of Ui when Ui(0) = 
u, it clearly holds that 

(7) ^sim(«l,«2) < tpa,nd(ui,U2) = 1pl(ui) + Tp 2 (u 2 ) - 1por{ui,U 2 ). 

As the ruin probabilities ip or , ip an d and ip s i m do not change under a scal- 
ing of (Ui,U 2 ), we restrict ourselves in the sequel to the respective ruin 
probabilities T/V^and and Vsim of the scaled process (X±,X 2 ) given by 
Xi(t) := Ui(t)/5i = Xi + Pit — S(t) with Xi = Ui/5i, pi = Ci/5i. This puts into 
evidence the fact that all the randomness in our model acts in one direction; 
in the future, we call this model informally a "two-dimensional degenerate 
model." 

Geometrical considerations. The solutions of the "degenerate two-di- 
mensional" ruin problems ip or , ?^ s i m and ^ an d strongly depend on the relative 
position of the vector of premium rates p = (j>i,p 2 ) with respect to the pro- 
portions vector (1,1). A key observation is that the ruin times r or and r s i m 
are also equal to 

(8) r 6 (xi, x 2 ) = mf{£ > : Sit) > b(t)}, 

where b = b m \ a = min{£i, £ 2 } and b = 6 max = max{4,^}, respectively, with 

(9) £i(t) = x, + p^, i = l,2,t>0; 

see Figure 1. The two-dimensional problems may thus also be viewed as a 
one-dimensional crossing problem over a (piecewise) linear barrier. 
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Our exact results follow directly from this geometrical observation, which 
essentially breaks the problem in two pieces: ruin of one of the coordinates 
before the deterministic time 

(10) T = T(x 1 ,x 2 ) = ^^± 

Pi ~P2 

of entering the lower cone x 2 < Xi, or ruin of the other coordinate subse- 
quently. 

If the initial reserves satisfy x 2 < x\, the two lines do not intersect. It 
follows therefore that the barriers b m i n , 6 max are actually linear 

b m in(t) =x 2 + p 2 t, b max (t) =xi +pit. 

Thus, the "or" and "sim" ruin always happen for the second and first com- 
pany, respectively, and 

1por(xi,X 2 ) =1p 2 (x 2 ), ■0sim(xi,X2) = 1pl(xi). 

This case in which explicit and asymptotic formulas for V'simiV'or and 
''/'and follow directly from the well-known one-dimensional ruin theory; see, 
for example, Rolski et al. (1999) or Asmussen (2000) — will not be discussed 
further. 
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Contents. This paper is devoted to obtaining exact results for the even- 
tual ruin probabilities in the other case, when the initial capitals satisfy 
x\ < X2- We obtain also sharp "Cramer /light tails" asymptotics, as the ini- 
tial capitals xi,x 2 tend to infinity along a ray (i.e., x\jx 2 is constant). 

We introduce notation and gather some necessary prerequisites from one 
dimensional ruin theory in Section 2. 

In Section 3, we obtain several exact decomposition formulas for the 
two-dimensional ruin probabilities Von Vsim arid ^and; in terms of one- 
dimensional ruin probabilities. In particular, we obtain for the compound 
Poisson case with exponential jumps an exact result in Corollary 1. 

In Section 4, we show in Theorem 3 for the general two-dimensional Levy 
case the natural result: 

■4>ot{x 1 ,x 2 ) ~^i(a;i) + ^2(^2), 

Vandal) x 2 ) = o(ipi(x 1 )+tp 2 (x 2 )), 

where we write f(x) = o(h(x)) (x — > 00) if lirn^oo f(x)/h(x) = and sim- 
ilarly f(x) ~ h(x) (x — ► 00) if lirn^oo f(x)/h(x) = 1. The corresponding 
renewal version is stated in Theorem 4. 

In Section 5, specializing to the case of the Cramer-Lundberg model, 
we sharpen the general result, obtaining two term asymptotic expansions — 
Theorems 5 and 6. We find different leading terms within subcones of the 
positive quadrant, as typical in such cases; see, for example, Borovkov and 
Mogulskii (2001) and Ignatyuk, Malyshev and Shcherbakov (1994). We also 
find in our specific "degenerate model" a correction term not present in 
previous works. 

The paper concludes with two explicit examples in Section 6, which pro- 
vides interesting illustrations of multidimensional sharp large deviations — 
see the Appendix for the relation to the existing first passage large deviations 
theory. In this context, it is worth clarifying that our particular results are 
considerably sharper than the general logarithmic asymptotics results ob- 
tained, for example, by Collamore (1996), Theorems 2.1 and 2.2, for the first 
passage times to convex open and open sets in M. d , respectively. Further, 
our proofs do not appeal at all to the multidimensional large deviations 
theory. Instead, exploiting the special structure of our problem, we make 
use of a one-dimensional asymptotic limit result; see Theorem 2 in Section 
2.3 — which is a consequence of the sharp approximation of finite-time ruin 
probabilities obtained by Arfwedson (1955) and Hoglund (1990). 

2. Preliminaries: One dimensional theory. Let Z be a general Levy pro- 
cess, that is, a process with stationary and independent increments that is 
continuous in probability and starts at 0, defined on some probability space 
(£l,J r , P) and let E denote the expectation w.r.t. P. To avoid trivialities, 
we exclude the case that Z has monotone paths. 
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We will restrict ourselves to Levy processes Z that admit negative expo- 
nential moments, that is, E[e~ uZ ^] < oo for some v > 0. For such processes, 
Z, we consider the cumulant exponent: 

k{6) :=t- l \og{E[e ez ®]), 

which is well defined on some maximal domain 

e = {9eR:K(9) <oo} 

whose interior will be denoted by 0° := (9, 9) , where 9 = inf and 9 = sup 0. 
The map 9 \— > k(9) restricted to the interval (9,9) is a convex differentiable 
function. By k'(9) [«+(#)], we denote the [left-] derivative of k at 9, respec- 
tively. 

In the particular case that Z is equal to the classical Cramer-Lundberg 
process (i.e., a positive drift added to a spectrally negative compound Pois- 
son process), we note that 9 = oo and k'(9) — > oo if 9 — > oo. 

2.1. One- dimensional Cramer type sharp asymptotics. Let r(x) denote 
the first passage time of the level — x by Z 

t(x) = inf{i > : x + Z(t) < 0} = M{t > : X(t) < 0}, 

where X is the translation of Z by x, X(t) = x + ^(i), and set ip{x) = 
P(t(x) < oo). Let us assume the Cramer assumption that there exists a 
7 > 0, such that «(— 7) = 0. 

Under this assumption, Cramer has shown that when Z is a Cramer- 
Lundberg process, then 

(11) lim e 7a Wx) = C, 

x — >oo 

where C is given explicitly by [see, e.g., Feller (1971), Chapter XII. 5] 

(12) C = -k'(0)/k' + (- 1 ), 

which is strictly positive if and only if /«+(— 7) > —00. For a general Levy 
process satisfying the Cramer assumption, Bertoin and Doney (1994) proved 
that (11) remains valid for some constant C > that can be expressed 
in terms of the law of the ladder process and that is positive precisely if 
E[\Z(l)\]<oo. 

2.2. The Arfwedson-Hoglund theorem. We recall now the asymptotics 
of the finite time ruin probability ijj(x,t) = P{t(x) < t) for the Cramer- 
Lundberg process, first obtained by Arfwedson (1955) via the saddle-point 
method. Later, Hoglund (1990) noted similar results for the probability of 
ruin after time t 



(13) w(t,x) =P(t < t(x) < 00) = if)(x) - ip(t,x). 
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Our formulation below, based on Hoglund (1990), Corollary 2.3 and Asmussen 
(2000), uses the exponential family of measures {P^} defined for all c G 6 = 
{9 : k(9) < oo} by the Radon-Nikodym derivative A( c ) 



dp( c ) 

(14) 



= A (c) (t) := exp(c(X t - X ) - n(c)t) = exp(cZ t - K(c)t) 

Tt 



dP 

with p(°) = P and the corresponding shifted cumulant exponent given by 

(15) k ( - c \9):=k(9 + c)-k{c). 

Further concepts, familiar from large deviations theory, that will be used 
are the convex conjugate k* of the cumulant exponent k, defined by 

(16) «*(«) = sup[«/3-f6(/3)], 



and the reparametrization of the exponential family by the corresponding 
set of means of Z. More precisely, it holds that to any —v £ (v,v) with 

v = ]hnK r (9) and v = lira k' (9) , 

019 e\e 

is associated a unique shift 9 E 6° = (6, 9) = (9, oo) such that 

k'(9 v ) = -v. 

Further, for any < v < —v, the conjugate shift 9' v € is defined via 

k(0 v ) = k(9' v ), 9 v <9' v 

and note that v > implies 9 V < 9' v . (Note that we parameterized by — v 
since the possible means leading from x > to are necessarily negative.) 

Remark 1. From now on, all quantities related to the shifted measure 
p( c ) (like, e.g., the cumulant exponent, the adjustment coefficient, etc.) will 
be indicated by a superscript (c) added to their P-counterparts. Observe 
that, in view of (15), the convex conjugate under P^ and P are related 
to each other by k*^(v) = k*(v) + k(c) — cv. Similarly, it follows that the 

shift 9^ and its conjugate 9^' are related to 9 V ,9' V by 9^ = 9 V — c and 

9^' = 9' v — c, respectively. 

Theorem 1. Assume that either k+(0) < or that the Cramer assump- 
tion holds and write ( = —mm{9:n(9) = 0}. If < v < —v and x,t — ► oo 
such that x/t = v; it holds that 

(17) Mxt)~{ Ce ~ X ^ ^/x/^<-<(-C), 
(i7j \\D(v)\t-y 2 e-^(-v\ if x /t>-K> + {-Q, 



(18) w(x,t) 



| J D(u)|r 1 /2 e -^*(-^) ) if x /t < - K ' + (-Q > 
Ce~<, l f x / t> - K > + (-(), 
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with C = 1 if C = and C = -k'(O) /k' + (-() if(>0 and 

(19) w-^-tssto -cW =w 

where (17) is to be understood as \im. x ^ OQjX= t v e x ^ij){x,t) = 0, i/ «'(—£) = 
— oo and C > 0. 

2.3. The asymptotic limit laws of the process before/after ruin. The fol- 
lowing result shows existence of and identifies the limit laws of the Cramer- 
Lundberg process X(t) conditioned on t < r and on t > r. It will be used to 
establish Propositions 2 and 3. 

Theorem 2. Assume that 9<0 and re'(0) < 0. 

(i) If0<v<-n'(0), then 

(20) % = lim P(X(t) G -lr(ar) > i), 

a;,t— >oo,a;=tD 

m £/ie sense of weak convergence, where 

(21) V v (dy) =c(v)- 1 [e-°°v - e^]l i0>oo) (y) dy 

with c(v) the function appearing in (19). 

(ii) If -re'(O) < v < -v, then 

(22) = lim P(X(t) e -|r(x) < t) 

as,t— >oo,£=ti> 

in i/ie sense of weak convergence, where 

(23) ^(dy) = | c ( w )r 1 [ e "^l ( o, 00 )(y) + e -^l ( „ 00i0) (y)] ( iy. 

Proof, (i) First, we verify that the measure ty v is a probability measure. 
Indeed, it is not hard to verify that ^f v is a measure (since 0' v > V > 0) that 
integrates to one. Further, it is easily checked that the mgf M v of is 
given by 

In view of the continuity theorem, the weak convergence in (20) follows once 
we show that the mgfs M x j(c) of the measures P(X{t) 6 • \t{x) > t) converge 
pointwise to M v (c) as x,t-> oo,x/t = v, for c in some neighborhood of the 
origin. 
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Since 9 < 0, it holds that for all c in some neighborhood of the origin 
M Xyt (c) is finite and E[e cX ^l {t<T{x)} ] = e ^+«(c)tp(c) ( t < r ( x ))_ It also holds 

that P^(t < t(x)) = P(°)(t < t(x) < oo) [since k( c )'(0+) < 0]. Therefore, 

M xA c) = EV*M\t < r(x)) = ^C^V^ 

P(t < t(x) < oo) 

P(t < t(x) < oo) 

Invoking Theorem 1 and (15) for the numerator and denominator, it follows 
by taking the limit of x, t — > oo, x = tv that 



lim M x>t (c)=D^(v)/D(v) 



}(c)'_/,(c) ,/,/. 



In view of Remark 1 , the latter is equal to M v (c) . 
The proof of (ii) is similar and omitted. □ 

2.4. Law of large numbers for the ruin time. We include now for refer- 
ence a result concerning the behavior of the time of ruin of a general Levy 
process Z for large initial reserves. 

Lemma 1. Suppose that E[\Z(1)\] < oo and E[Z(1)] < 0. Then, as x — > 
oo; 

(i) r(x)/s-> -EiZil)]- 1 P-a.s. and 

(ii ) E[t{x)]/x-+-E[Z(1)]-\ 

Proof, (ii) If i?[Z(l)] = 0, then the Levy process Z oscillates and the 
identity follows since then E[t(x)] = +oo for every x (see, e.g., Bertoin 
(1996), Chapter VI, Proposition 17(iii)). Suppose now that — oo < £ ! [2'(1)] < 
(so that Z drifts to — oo) and first exclude the case that Z is a compound 
Poisson process. Denoting by L~ l {t) = inf{u > 0:L(u) > t}, the inverse of 
the local time L of Z and T(x) = inf{t > : H (t) > x} the first passage time 
of the ladder height process H{t) = Z(L -1 (t)) of Z it is easily verified that 
t(x) = L~ l (T(x)). The pair (L _1 ,i?) forms a two-dimensional Levy process 
and we denote its bivariate Laplace exponent by k. The Laplace transform 
of -E[r(x)] can then be expressed as follows: 

i r „ r , Nn , 5ik(0,0 + ) 

where di denotes the partial derivative with respect to ith variable (see, 
e.g., Bertoin (1996), Chapter VI, Proposition 17). As £(0,0) =0 and k(0, •) 
is right-differentiable in zero, it follows in view of a Tauberian theorem that 

E[t(x)] ~x<9ik(0,0 + )/<9 2 £(0,0 + ) as x^oo. 
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The strong law of large numbers implies that the product H(t)/L 1 (t) = 
x [t/L^it)} converges to 

(24) E\Z(l))=E[Z(L-\l)))E[L-\l)r 1 

(the corresponding result for random walks is known as the famous Wald 
identity). Since <9i£(0+,0) = E[L~ l {l)} and d^,^) = E[Z(L~ l {!))}, the 
claim follows. The case of a compound Poisson process follows by adding a 
small drift. 

(i) The strong law of large numbers implies that, P-a.s., 
t{x)/x = L-\T{x))/T{x) ■ T(x)/x - E[L~\l)}/ E[H{1)\ = E[Z(1)], 
as x — > oo, where we used the Wald identity (24). □ 

3. The exact ultimate ruin probability for the degenerate 2-d process. 

In this section, we consider the probability that a Levy process S starting 
at ever upcrosses a piecewise linear barrier b. To be specific, we consider 
the first passage time 77, of S over b, as in (8), where b is given by b = 6 m i n 
or b = 6 max with 

&min(i) = mm {x i + pit} , b max (t) = max{xj +pit} 

where X2 > x\ and pi > p 2 . As noted in the Introduction, P(r& < 00) with b = 
&min (resp. b = b max ) exactly coincides with the ruin probability ip or (xi, X2) 
[resp. ip s i m (xi,X2)] of the process (Xi,X 2 ) with 

(25) Xi(t) :=Xi + Pit- S(t), i = l,2. 
Denoting by 

(26) T = T(xi,x 2 ) = ^—^ 

Pi ~P2 

the time at which the lines 1 1— > x\ +p\t and 1 1— > x 2 +P2t cross, we see that, for 
example, for S to never cross b m \ n , it is required to stay below the barrier 
x i + Pit between the times and T and subsequently to stay below the 
barrier x 2 -\-p2t after time T. Since S is Markovian and x\ +p{T = x 2 +p 2 T, 
conditioning at time T yields 

coo 



poo 

ip OT (xi,x 2 ) = fi(ds,T\x 1 )F 2 (s) 
Jo 



10 

where ^ or = 1 - tp or , Fi(s) = P(s + S(t) <Xi+pi(t + T) Mt > 0) and 

(27) fi{ds,T\x):=P{S(t)<x + Pl t, Vt G [0, T],S(T) € ds) 

is the density of S(T) of the paths at time T that "survived" the upper 
barrier x + Pit. Reformulating this result in terms of the two coordinates Xi 
of X = (Xi,X 2 ) and the coordinate- wise densities of the surviving paths 

(28) ^i(dz,T\xi) = P{n(xi) > T,Xi(T) G dz) 
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we arrive thus at the following result, which relates the ruin probabilities of 
the two dimensional process X to those of its coordinates X\ , X2 ■ 

Proposition 1. Let X = (Xi,X 2 ) be the two-dimensional Levy process 
with Xi given in (25) and suppose that x 2 > x\ and p 2 <Pi- 

(a) The ruin probabilities Vwi j V'sim an d "0or ar e given by 
(29) Vsim(xi, x 2 ) = P(r 2 (x 2 ) < T) +p(t 2 (x 2 ) > T, inf X 1 (s) < 

\ s>T 



(30) ^or(xi,x 2 ) = P(n{xi) < T) +P[t 1 {xi) > T,MX 2 (s) < 

(31) VwOci,^) = P(T < ri(xi) < 00) + P{t\(x\) < T,t 2 (x 2 ) < 00). 

(b) The survival probabilities ip OI = 1 — tp or and tp s { m = 1 — tpsim are given 

by 

ip sim (x 1 ,x 2 ) = P\r 2 (x 2 ) > T, inf X^s) > O^j = tp^dz^x^ip^z), 

■0 or (xi,x 2 ) =p(n(xi) > T, inf X 2 (s) > O^j = ip^dz^x^^z), 

where T is given in (26), V^(^-z, T|xj) in (28) and tpi(z) = P{ti{z) = 00) are 
perpetual one- dimensional survival probabilities. 

Proof. By definition of ^sim, it holds that 

^ sim (xi,2 2 ) = P(max{Xi(t),X 2 (t)} > for all t > 0). 

Next, we note that, if x 2 > x±, it holds that the maximum 

max{Xi(t),X 2 (t)} = max{xi - x 2 + (pi -p 2 )t,0} + X 2 (t) 

is equal to X 2 (t) for t <T and to Xi(t) for t > T, where T was defined in 
(26). Applying subsequently the Markov property of X\ at time T shows 
that 

VwO^l^) = P{X 2 {t) > for t < T,Xi(t) > for t > T) 

(32) 



/ P(X 2 (T)edz,r(^)>T)P{T 1 (z)= 00). 
JO 



The identity in (29) follows by taking the complement of (32). The proof of 
ip or is similar and omitted. Finally, write 

^and(2T,£2) = P(T < Tl(.Xl) < 00,T 2 (x 2 ) < Oo) + P(ti(xi) < T,T 2 {x 2 ) < Oo). 

Equation (31) follows then by checking that {t\(xi) > T} and {inf s> T Xl(s) < 
0} respectively imply that {r 2 (x 2 ) > T} and {ini s> T X 2 (s) < 0}. □ 
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In the special case that S is a compound Poisson process with exponential 
claims <7j with parameter fj,, we have exponential ultimate ruin probabilities 

(33) ^ i ( x ) = C l e-^ x , 

where Cj = -^7 and 7^ = fi — X/pi. Similarly, if S is a spectrally negative Levy 
process, the Markov property and the absence of positive jumps imply the 
multiplicativity property P(ji(x + y) < 00) = P(ti{x) < oo)P(rj(?/) < 00). 
Thus, P(Ti(x) < 00) must be an exponential function (33) and the constant 
Ci equals 1. For these two cases, equations (29)-(31) can be developed fur- 
ther by employing the technique of change of measure. 

3.1. The case of exponential ultimate ruin probabilities. Consider now 
the relation (30) in the case of exponential ultimate ruin probabilities that 
is when ipi(x) = de~^ x for d,^ > (i = 1,2). Note that 

(34) Vor(zi,x 2 ) = P(ti(zi) < T) + C 2 E[e~^ T h {Tl{xi)>T} }, 

where C2 = 1 in the case that S is a spectrally negative Levy process. Let 
Ki be the cumulant exponent of Xi(t) — Xi. By a change of measure (14) and 
using that —72^1 + 72 )^ = —72^2, we find that the second term in (34) 
is equal to 

C 2 e-^ xl+K1 ^ T E[A^\T)l {Tl{xi)>T} ] = C 2 e-^ X2 P(-^(T 1 (x 1 ) > T). 

The probabilities Vsim can be treated using similar arguments, and Vand 
is obtained from the "complementarity equation" (7). In conclusion, the 
original two-dimensional ruin problems V ; or/''/'sim/^ ; and ar e reduced to one- 
dimensional finite time ruin problems f/y (x, t) = p( c \ri(x) < t), as follows. 

Corollary 1. Suppose S is a spectrally negative Levy process, or a 
compound Poisson process with exponential (positive) jumps and let Xi de- 
fined by (25). If X2 > x\, it holds that 

Aim(xi,X 2 ) = 1p 2 (x 2 , T) + Vl (£l)V4 11] {X2,T), 
^or (2:i,a:2) =iJl(x 1 ,T) + $ 2 {x 2 )^ x " /2 \x 1 ,T), 
^and^l,^) = Wl{x X ,T) + Tp 2 {x 2 )Tp[~ l2) (Xi,T), 

where i^f 1 (x,t) = 1 — ijJf\x,t) and 

(35) wi (x,t) = P(t < n (x) < 00) = ipi (x) -tpi(x,t). 
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Proof. Let us establish the last statement. By inserting the expression 
for ip or in (7), it follows that 

VWdf^l,^) = "01 (Zl) +^2(^2) -1por(xi,X 2 ) 

= Mxi)-Mxi,T)+Mx2)(l-^[~ l2 \xi,T)). □ 

This decomposition result (and its generalization) will provide the key for 
obtaining the two terms asymptotic expansions in Propositions 2, 3 below. 

4. General two-dimensional Cramer asymptotics. We consider now the 
asymptotics of the ruin probabilities ip s im, and ipa.nd when the initial re- 
serves tend to infinity along a ray, for a general two-dimensional Levy process 
X = (Xi, X2) starting from x = {x\,X2). To avoid degeneracies, we exclude 
throughout the cases that X± or X 2 have monotone paths, or that the ratio 
\X\ — x{\j\Xi — X2] is constant. The law of the process X is determined 
by its joint cumulant exponent k(0i,0 2 ) = log£[e ei(Xl(1) ~ Xl)+e2(X2(1)_a:2) ] 
which is well defined on its domain S = {0 G M? : k{9) < 00}, whose interior 
is denoted by 3°. For every G S°, the gradient V«(0) = (din(6), 82^(6)) is 
well defined. Other subsets of S playing a role in our setting are: the Cramer 
set C, its interior C° and its boundary dC := C\C° where 

C = {(0i,0 2 )GH:«(0 1 ,0 2 )<O}. 

In view of the convexity of k, it follows that the set C is convex and that for 
fixed 0' G dC n E°, it holds that 

(36) [0 - 0'] • Vk(0') < for all G C, 

where • denotes the inner-product. 

Associated to any c G S is a measure defined twist of P by the 
martingale exp(ci(Xi(i) — x±) + c 2 (X 2 (t) — rr 2 ) — k(ci,c 2 )£). 

We assume throughout that besides the origin, the Cramer set intersects 
the axes in two more points: 7W = (—71,0), 7^ = (0,-72) G c?C,7i,7 2 > 0, 
so that 

(37) «(- 7 i,0) = k(0,- 72 ) = 0. 
We shall also assume that 

(38) (- 71 ,0),(0,- 72 )GH°. 



When (37), (38) hold, we will say that the Cramer assumptions hold true 
for X = (X 1 ,X 2 ). 
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Example 1. If S(t) is a Levy process and Xi(t) = Xi + pit — S(t) (with 
the cumulant generating functions Kj), then the joint cumulant generating 
function k, of (X\,X 2 ) is related to K\ and k 2 by 

n{e u e 2 ) = + e 2 ) - e 2 (pi - p 2 ) = k 2 (0i + 62) + 0i(pi - &) 
=p-9 + K S (0i + e 2 ), 

where ks is the cumulant exponent of S. It is easy to check that the de- 
generate two-dimensional Levy process X = (X\,X 2 ) satisfies the Cramer- 
conditions iff its coordinates do, that is, if there exist constants ji > 0, in the 
interior of the domains of the cumulant exponents Ki of Xi(t) — Xi (i = 1, 2), 
such that 

(39) Ki(-7i)=0. 

The following result yields the asymptotics of ip or and an order estimate 
of Vwa for general two-dimensional Levy processes. 

Theorem 3. Suppose that the Cramer assumptions (37) and (38) hold, 
and let a > 0. Then as K — ► oo, 

(40) ^ or (aK, K) ~ C 2 e~^ K + C x e~ lxaK ', 

(41) ^ nd (aK,K) = o(C 2 e-^ K + C ie -^ aK ), 

where Ci > 0,i = 1,2, are i/te asymptotic constants corresponding to Xi. 

The proof of Theorem 3 is based on the following estimates. 

Lemma 2. The following hold true: 

(i) max{ipi(xi),ip 2 (x 2 )} < ^ r(xi,x 2 ) < tpi(xt) +ip 2 (x 2 ); 

(ii) ipor(xi, x 2 ) =tpi< 2 (x 1 , x 2 ) + tp 2 <i(x 1 ,x 2 ) - ipi= 2 (xi,x 2 ), where 

tpi<j(xi,X 2 ) \=P{Ti(xi) < Tj(Xj),Ti(Xi) < oo) 

and 

^1=2(^1,^2) :=-P(ri(xi) =t 2 (x 2 ) < 00). 

Proof. The estimates follow in view of the observations that 

2 

{Ti(xi) < 00} C {t ot (x 1 ,x 2 ) < 00} C {jinixi) < 00} for i = 1,2, 

i=l 

and 

{r or (xi,x 2 ) < 00} = A x U A 2 \[Ai n A 2 ], 
where A; = {n{xi) < 00, r^Xi) < T 3 _j(x 3 _j)}. □ 
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Lemma 3. Suppose that (37) and (38) hold, and write 7 0j/ g := flaji + 



(42) Vi<2(aK,K)~Cie-^, ^2<i (aK, K) ~ C 2 e~^ K . 

(ii) For a>0 and any (3 G (0,1), ip sim (aK , K) = o(e~^ K ) (K->oo). 



This lemma (established below) implies immediately Theorem 3. 

PROOF of Theorem 3. First note that, in view of the Cramer-Lundberg 
asymptotics (11) and equation (7), the asymptotics in (40) imply the esti- 
mate in (41). The rest of the proof is therefore devoted to establishing (40). 

In view of (11) and Lemma 2(i), it follows that, if 71a > 72 [resp. 71a < 72], 
the lower bound and upper bound in Lemma 2(i) are of the same order of 
magnitude, C 2 e~ 12K [resp. C\e~' yiaK ], as K — > 00. Thus, (40) is valid if 
71a 7^72. 

Next we turn to the case 71a = 72. Since ■01=2 is dominated by ip s i m 
and 7 a/ g = 72 = aj± if 071 = 72, it follows, by invoking Lemma 3(h), that 
ipi =2 (aK, K) = o(e~" l2K ) = o(e~" liaK ) as K — > 00. In view of Lemma 2(h) 
and Lemma 3(i), it therefore follows that (40) is also valid if 071 =72. □ 

Proof of Lemma 3. (i) The asymptotics of ipi< 2 follow once we have 
shown that as K — > 00 it holds that 



(43) e^ aK ^< 2 (aK, K) = E^\e-^ x ^l {ri ^ Tl<oo] ) - C u 



where T\ = T\(aK) and t 2 = t 2 (K). To prove this claim, we compare the 
asymptotic behavior of t\ and t 2 as K — > 00, adapting the argument devel- 
oped in Glasserman and Wang (1997) (Proposition 2) for random walk. If 
£(-7i,o)[x 2 (l)] > X2 , then p(-^' )(r 2 = 00) -> 1 as K -> 00 and, invoking 
(11), the claim (43) follows. If E^ 1 ^ [X 2 (l)\ < x 2 , it follows in view of 
Lemma l(i) that as K — >oo, p( _7l '°)-a.s., 



where we used that 8^(6) = E^[Xi(l) - asj for 9 G 3°, i = l, 2. Applying 
(36) with 6 = (0, —72) and 6' = (—71,0), we see that the right-hand side of 
(44) is bounded above by 071/72, which is equal to one if 72 = 071. There- 
fore, t 2 (K) dominates T\(aK) for all K large enough and (43) follows as a 
consequence of the Cramer-Lundberg asymptotics (11). The asymptotics of 
-02 <i can be treated similarly. 



(ii) Choose /3 € (0,1) and write 7(^ = ^(71,0) + (1 -£)((), 72) = (71 (/?), 
J 2 (P))- The key step is to verify that the segment {7^), /3 G (0,1)} is not 




(44) 



Ti(aK) 
r 2 {K) 



= a 



n(aK) K 
aK t 2 (K) 



ad 2 t(-7i>0) 
9i«(-7i,0) 
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part of the boundary dC. Indeed, since the function (3 t— > f{(3) with f(f3) = 
exp(/e(- 7 (/?))) = f?[exp{/3 7 i(Xi(l) - x x ) + (1 - /3)7 2 (X 2 (1) - s 2 )}] is strictly 
convex with /(0) = /(l) = 1, it follows that /(/?) < 1 for f3 G (0,1). (The 
strict convexity is a consequence of the facts that /"(/?) = £'[(7i(Xi(l) — 
a*) - 72(^2(1) - x a )) 2 e^CXi(i)-*i)+(iH3baCX2(i)-*a)] > and we excluded 
the case [X\ — X\\/\X 2 — x 2 ] = const.) 

Therefore, there exists a —7* = —(71,72) £ C° such that 7* > 7i(/3) (i = 
1,2). By changing the measure, we see that tp s i m (aK,K) is equal to 

p -(7?o+7|)K R (-7*)r 7|Jf 2 (r sim )+75 , Xi(r sim )+K(-7j , ,-7|)r sim1 1 
c ^ L c J -{T sim <oo}j! 

where r s i m = r B i in (aiif, if). Since Xj(r s i m ) < and «(— 7*, —72) < 0, this ex- 
pectation is bounded above by 1, and as 07 j" + 7^ > 7 a ,/3, it thus follows that 
ip sim (aK,K) =o{e~^ K ) as K^oo. □' 

The following result concerns the asymptotics in the upper cone {x± < 
X2}, in the case of the Sparre-Andersen model. 

Theorem 4. Let S be a compound renewal process as in (2) and let 
a < 1. Assume there exist 7i > and an e > s-uc/t £/iai I?[e _7iP ^]i?[e 7iCT ] = 1 
and E[e^ i+t ^ a ] < 00. Then it holds that, as K — > 00, 

(45) 1><x(aK, K) ~ f^e" 72 ^ + de" 71 ^, 

(46) Vand(aK, K) = o(C 2 e- 72X + C x e~^ aK ). 

Proof. In view of the key observation that the ruin probabilities s/w/V'and 
and Vsim do not change if we replace X = (X\,X2) by a two-dimensional 
compound Poisson process with unit jump rate and jump sizes distributed 
as (<7 n — piCn,o~n — P2Cn) , the statement follows by invoking Theorem 3. □ 

In Section 5 below, we will sharpen Theorem 3, in the degenerate case. 
Before that, we introduce a partition of the quadrant in cones, which turn 
out to describe the different asymptotic regimes of the ruin probabilities as 
the initial reserves (x\,X2) tend to infinity along a ray. 

4.1. The asymptotic cones. We introduce now two cones T>i (i = 1,2) 
within the quadrant = (0,oo) 2 , situated between the and the 

directions of the expected drift evaluated at the adjustment tilts vW := 
Vft( 7 ^),i = 1,2. Heuristically, these cones can be described as the "asymp- 
totic boundary cones" of the "sim" ruin, that is, the cones where asymp- 
totically this event happens dominantly by straight paths running to the 
half-lines {x : xi = 0, x%_i < 0, i = 1, 2}. 

Letting denote the "slopes ^1" of the vector Vk( 7 W) (i = 1,2), the 
cones are given by 

T>i = {(xi,x 2 ) Sl^:xi> x 2 si}, 



EXIT PROBLEM OF A TWO-DIMENSIONAL RISK PROCESS 



17 



(47) 

£>2 = {(x\,x 2 ) G M. 2 + :xi < x 2 s 2 }. 
Note that in the degenerate case, Si become 

(48) *i = 4£ 3 4. s 2 -( K '^ 



4(-7i)' " V4(-72) 

and we show in Lemma 4 that s 2 < s±, implying that the cones T>\ and T> 2 
are disjoint in this case. Let now 

V = Rl\\V 1 UV 2 } 

denote the open cone lying between T>\ and T> 2 (where T>i is the closure of 
Vi). 

We will show in Theorem 5 that within T>i, ip s i m (xi,x 2 ) is asymptotically 
equivalent to tpi(xi),i = 1, 2, respectively, and that a different regime holds 
within Dq , which is characterized by "radial dependence" on the slope a = 

We will also show in Theorem 6 that in the case when V 2 is void [which 
is characterized by diKi(j^) > 0], a special type asymptotic regime holds 
for Vand) within a new "secondary cone" T> 2 situated between the x 2 axis 
and the direction v^ 3 ) := Vk^ 3 )), where is defined as the leftmost 
intersection of the Cramer set with the line 9 2 = — r y 2 . In the degenerate 
case, we have 

£>2 = {(xi,x 2 ) £ R 2 + :x\ < x 2 s 3 } 

where 

S3 = 4(-73)/4(-73) 

with 73 the largest root of k±(—s) = 72). As stated in Lemma 4 below, 
if k'i (—72) > 0, it follows that T> 2 7^ = T> 2 and otherwise the cones T> 2 and 
T> 2 coincide. This partition of the positive quadrant into cones 

Vt,V 2 and Po:=K+\fPiUP 2 ] 
is illustrated in Figure 2. 

5. Sharp asymptotics for degenerate risk-processes. We restrict now 
ourselves to a two-dimensional Levy process (X\,X 2 ) with Xi(t) = Xi+pit — 
S{t), i = 1,2, where S is a compound Poisson process with positive jumps. 
Throughout this section, we assume that 

there exist 7« > 0, i = 1, 2, such that «»(— 7i) = 0, 

where Ki is the cumulant exponent of Xi — Xj. Note that whenever 71 exists, 
72 exists as well, and 71 > 72 (by the convexity of Ki). 
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FlG. 2. Pictured are the positive quadrant, divided in the three cones D\ (shaded with 
vertical lines), T>2 or D2 (shaded with horizontal lines) and Do or Do (white), as well as 
the lines Xi = X2 (dashed) and 71X1 = 722:2 (dotted). 

5.1. A characterization of the asymptotic cones. We start with gathering 
some properties of the asymptotic cones, in the degenerate setting. 

Lemma 4. The following hold true: 



(i) The cones T>i,i = 0, 1,2 are disjoint and T>q,T>i 7= 0. 

(ii) V 2 = V 2 + iff <(-72) < and V 2 = + V 2 iff < (-72) > 0. 

(iii) VxrzU:={(xi,x 2 )eRl:x2j2<x l -fx} andV 2 cR 2 + \U. 

Proof. Writing 



it follows that si < 1, since k' 2 {— 71) < 0, and that s 2 < s±, since 71 > 72 and, 
by the strict convexity of k 2 , k 2 is strictly increasing on its domain. Next, 
in view of the definitions of s 2 and S3, it follows that S3 = [resp. S2 = 0] 
iff k'i(— 72) = [resp. 72) > 0]. Subsequently, we note that on the ray 
x i/ x 2 = 72/71 it holds that 

x 2 P1-P2 = «2(~7i) - «d(-7i) = «2(~7i) -^(-72) 

T(x 1 ,x 2 ) 1- 72/71 71 -72 71 -72 

The strict convexity of k 2 thus implies that along the ray x\jx 2 =72/71 it 
holds that — n' 2 {— 72) < x 2 /T(x\,x 2 ) < — K 2 (— 71). It is a matter of algebra 
to verify that these inequalities are equivalent to S2 < 72/71 < si (see also 
Lemma 5 below). The assertions (i), (ii) and (iii) follow then in view of the 
definitions of T>i, i = 0, 1, 2, and T> 2 . □ 



k'i{s) = k! 2 {s) +pi-p 2 
4(s) k' 2 (s) 



1 + 



Pi ~P2 
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A key point in the analysis of the asymptotics of the degenerate risk 
processes is an equivalent description of the cones T>i in terms of comparisons 
with the time T = T(x±,x 2 ) defined in (26), which will enable us to translate 
the asymptotics in two-dimensional space into the "space-time" asymptotics 
of Arfwedson (1955) and Hdglund (1990). 

Lemma 5. Writing Ti = Xi/[— 7i)] and T{ = Xi/[— k[(— 73-i)], i = 
1,2, and a(x) = ^ [where x = (x±,x 2 ) G the following hold true: 

V 2 = {x£R 2 + : a(x) < s 2 } = {ieR+:fi<T} = {iGl 2 + :T 2 < T}, 
D = {xei+:Ti<T<T 2 }, 

V 1 = {xeM 2 + :s 1 < a(x)} = {x G M. 2 + : T < f 2 } = {x G M. 2 + : T < T x }. 

Proof. The first two equalities for T>\ and T> 2 are just the definitions, 
given in (47). Next, we note that from the definition of T it is easy to check 
that 

(49) xi/x 2 = a 44> x 2 /T(x 1 ,x 2 ) = (pi-p 2 )/(l-a)=v a . 

In particular, inserting a = si [defined in (48)] and using k'i(s) = pi — p 2 + 
k' 2 (s) it follows that s± = x\jx 2 iff T = x±/[— 71)] iff T = x 2 /[— k' 2 (— 71)]. 
The two representations for T>\ now directly follow. The identities for T> 2 are 
proved similarly. Finally, the equalities for T>q follow by intersecting those 
for the complements T> 2 and T> 1 . □ 

5.2. Asymptotics. The leading term asymptotics of the two-dimensional 
ruin probabilities will be expressed in terms of the usual "one dimensional 
large deviations cast": the adjustment coefficients 7i > of Xj satisfying 
7i) = 0, and 7(a) given for < a < 1 by 

7 (a) = K* 2 {-v a )/v a where v a := (pi -P2)/{1 - a). 

Below we consider asymptotics along the rays (aK, K) in the plane with 
a <a, where 

a=l + (p 1 -p 2 )/v 

with, as before, 

(50) e = mi{6:K 2 (9) < 00} and v = lim k' 2 (6). 

919 

Note that in terms of the "space-time velocities" the restriction a <a reads 
as v = v a < —v. 
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Theorem 5. Assume that 9 < —71. If a < a, it holds as K — > 00, 



c ie -^ aK , 

^j sim (aK,K) ~ { (D*(v a ) + D' 2 (v a ))K- l / 2 e-^ K , 
C 2 e~^ K , 

where, for i = 1,2, Ci = — k^(0)/k-(— 7») and 
(51) = ^ 



if(aK,K)eV 1 , 

if(aK,K)eV , 
if(aK,K)eV 2 , 



w 



(52) /#(«,) 



1 



1 



+ 



«3-i(o) 4-i(o) 



where 6 W < 9w satisfy k' 2 (9 w ) = —w and Ki(6 w ) = Ki(9w )■ 

Next, we turn to the asymptotics of the ruin probability VWb which are 
formulated in terms of the (sometimes different) partition of the positive 

quadrant into V 1 , V 2 and V = U V 2 \. 

Theorem 6. If a<a, then it holds that, as K — > 00, 

(C ie -^ aK , if(aK,K)£V 1 , 

^ nd (aK,K) ~ J {D'^va) - D*{v a ))K- l l 2 e-^ K , if (aK,K) eV , 
{ C 2 e-^ :i+ ^- a ^ K , if {aK, K) G V 2 , 

where 

C 2 = -C 2K / 1 (- 72 )K(-73) 
and D[ and Df are respectively given by (51) and (52). 

Note. These results imply that if (aK, K) is contained in either T>\ or 
T>2, then tjj S i m (aK,K) and tp^diaK, K) are of the same order. 

5.3. Two terms asymptotic expansions in terms of one- dimensional shifted 
measures. The key for obtaining "two terms asymptotic expansions" for ip or 
and ipsi m (and then also leading term asymptotics) is given by the following 
decompositions that are generalizations of Corollary 1 to the current setting. 



Corollary 2. Let x 2 > x\. It holds that 
iJj t(xi,x 2 ) = ipi(x!,T) + C 2 (xi,T)e 



ihm(xi,x 2 ) = M*2,T) + d(x 2 , T)e~^ x ^ 2 n '(x 2 , T) 



72) (X!,T), 



EXIT PROBLEM OF A TWO-DIMENSIONAL RISK PROCESS 21 

where, for i = l,2, 

C i (x 3 - i ,T) = E^[h i (X 3 - i (T))\T3-i(x 3 - i )>T] 
with hi(x) = e' 1iX il)i(x),i = 1,2. 

Proof. In view of (30), the Markov property and a change of measure 
it follows that 

(53) ^ ot (x 1 ,x 2 ) = Mxi,T)+p(t 1 (x 1 )>T, inf X 2 (s) < o] 

\ T<s<oo / 

(54) = Vi (xi , T) + e-^ E(~^ [h 2 (X 1 (T))1 {ti (sbi)>t} ] 

using that A"i (T) = X 2 (T) . The proof of the decomposition of ^sim is similar 
and omitted. □ 

For V'and similar decompositions are derived in the following result. Recall 
that 73 is defined as the largest root of n\ (—72) = k\{—9) and set 7 := 73 — 72- 

Corollary 3. For x 2 > x\ it holds that 

(55) ^sad(xi,x 2 ) =wx(x 1 ,T) + P(ti(x 1 ) <T,t 2 {x 2 ) < oo) 
where w\ is given in (35) and, with n = t\{x\) and t 2 = t 2 (x 2 ), 

e^+^Pin <T,t 2 < oo) 

(56) =^-^-^)[e^ 1 ^)/ l 2(X 2 (T 1 ))l {T1 < T} ] 

= e^C 2 (x 2 ,T)4-^\x 2 ,T) + C 1 (x 1 ,T)4-' (3 \x u T)^ 
forC 2 (x 2 ,T) = E^-^[e^ x ^ T ^ 2 (X 2 (T))\T 2 < T] and 

C 1 (x 1 ,T)=E^[e^ x ^ T H 2 {X 1 {T))\T 1 < T\. 

Proof. Recall that (55) was derived in Proposition 1. It follows by def- 
inition of 73 that 72 — 7) = 72) or, equivalently, k[ 72 - ) (— 7) = 0. 
In view of this observation and the form of k(u,v), derived in Example 
1, it follows that k(— 7,— 72) =0. Changing measure with the martingale 
exp(— 7(A"i(t) — x\) — j 2 (X 2 (t) — x 2 )) and applying the strong Markov prop- 
erty at 77 yields the first equality in (56). 

The second equality follows by noting that 

P(n <t,t 2 < 00) 

= P(n <t,t 2 <t) + P( Tl < t,t < t 2 < 00) 
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= P{T2 <T)+ E[1 {ti < t<T2} P X2{t) (t 2 < co)] 

= E[l {T2 ^ 2 (X 2 {T))]+E[l {T ^ T} i, 2 {X 2 {T))\ 

= e-^ X2 E^[e^ X2iT ^ 2 (X 2 (T))\T 2 <T]tp!f 12 \x2,T) 

+ e -72x 2 -7xi x E^\e^ x ^ T ^ 2 {Xi{T))\T l <T)^[- ri \x l ,T), 

where in the second line we used the Markov property and that {t 2 < T} C 
{t\ < T} and in the last line we changed the measure and that X\ (T) = 
X 2 (T). □ 

We write 

f~g + h as x — > oo iff lim (/ — g)/h(x) = lim (/ — h)/g(x) = 1. 



Proposition 2. Assume that 8< -71. For any v < —v, it holds that, 
if xi, x 2 — > 00 such that x 2 /T(xi,x 2 ) = v, 

(57) (xi , x 2 ) « Vi (*i , T) + C 2 {v)e-^ x ^[~ l2 \x l , T) , 

(58) Aim(x u x 2 ) ^^ 2 (x 2 ,T) + Ci(i;)e- 7lXl 4" 7l) (x 2 ,r), 
where, fori = 1,2 and v / — k' 2 (— 71), — k' 2 (— 72), 

<W \ C3 _i(«,7i)- 1 [^(^) " #$ 3 ~°)], »/0 < « < -4(-7i), 
where 9 V < 0« satisfy k' 2 (9 v ) = —v, Ki(0 v ) = Ki{9y^), and 

/ \ u v — t/ v 

Ci(V,C) = 7TT 

with ij}* being the Laplace transform offy. 

Proposition 3. For any v < —v, it holds that, if x\,x 2 — > 00 such that 
x 2 /T(xi,x 2 ) = v, 

^ a nd{xi,x 2 ) 

(59) «Vi(zi)-V>i(zi,T) 

+ {C 2 ( W )e-^ X2 4" 72) (x2,T) + C 1 (t;)e-^ a;2 -^ 1 4" 73) (xi,T)}, 
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where 7 = 73 — 72 and for v / — k' 2 (— 73), 

0, ifO<v<-K' 2 (-j 3 ), 
|c 2 K72)|- 1 -^(#), ifv>-n' 2 (- l3 ), 

-C 2 4t = ^, ifO<v<-K> 2 (- 73 ), 

K l(-73j 

ci^^s)!" 1 • me { v ] ) - e- 1 ) ifv> -4(-7 3 ), 

with ifl being the Laplace transform of . 

Proof of Proposition 2. In view of Corollary 2, the proof for ip OI 
is complete once we show that C 2 (xi,T) converges to C 2 (v) if xi,x 2 — > 00 
such that X2/T = v. We distinguish between two cases. 

If x\ + k'i(— 72)^ > 0, then the strong law of large numbers implies that 
X ± (T) -> 00 p(-T 2 )-a.s. and that P^Xt^x-l) < T) tends to zero (see The- 
orem 1). Since h 2 (y) — > C2 as y — ► 00 (by the Cramer-Lundberg approxima- 
tion), we conclude that C 2 (xi, T) converges to C 2 (by bounded convergence). 

In the case that x\ + (— 72 )T < 0, we note that by virtue of Theorem 2(i) 
the distribution of Xi(T) conditioned on t\>T (under P^~^) converges to 
the measure *$> v in (23), with the shifts in (23) calculated using the cumulant 

exponent k[ 72 ^ and direction v' = x\/T(xi, x 2 ). Note that the application 
of Theorem 2 is justified since we have that 9 < —71 < —72. Indeed, observe 
at this point that k\ 72 ^'(#i 72 ^) = k'i(0 v ) = —v' and that, in view of Remark 
1, 9 ( f l2) = 9 V + 72 and ^ 1)( ~ 72) = ei 1] + 72 . Thus, C 2 (x 2 ,T) converges to 

{9 v + l2 ){9^+ l2 )^. v __«W t 



(60) / h 2 ( y )* v (dy)= My) " nr ^ [ e -**- e -»« 

The proof of the asymptotics of ip s i m is similar and omitted. □ 



Proof of Proposition 3. In view of Corollary 3 to finish the proof, 
we have to show convergence of the conditional expectations in the two 
different cases. 

In the first case when x\ + k'i(— 73 )T < 0, it follows by the law of large 
numbers that p(~ 73 )(ri < T) tends to 1. Also, taking note of Lemma 1 
and of the fact that in view of the definition of T, it holds that A"2(ti) = 
X\{T\) + (pi — P2)[T — n]), it follows that X 2 (n) — > 00 and h 2 (X 2 (ri)) — > C 2 , 
p(-73)_ a s _ Therefore, the bounded convergence theorem implies that the 
expectation in the first line of (56) converges to C 2 C (where C denotes the 
asymptotic constant for tpi under p(~T 2 )). 

In the opposite case that x\ + k'i{— 73 )T > invoking Theorem 2 as in 
Proposition 2 (which is in this case justified as 9 < —71 < —73 < —72) yields 
the form of C\{v) and C 2 {v) and in view of (31), the proof is complete. □ 
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5.4. Proofs of Theorems 5 and 6. In the proof, we use the following 
result: 

Lemma 6. (i) We have 7(a) > max{o7i, 72} for a ^ s±, S2 and 

7(0) = 071 + «2 ( _w o)/«a = 071 + Ki^'i-avJ/Va = K*(-av a )/v a . 

(ii) //k-(0+) > 0, it holds that ip*(9) = 9' 1 - «$(())/«< (0). 

Proof, (i) is a direct consequence of the definitions of 7(a), % and Re- 
mark 1. (ii) directly follows from Bertoin (1996), Theorem VII. 10. □ 

Proof of Theorems 5 and 6. In view of (58), the proof of Theorem 
5 consists in identifying the leading order term by applying the Arfwedson- 
Hoglund's Theorem 1 to the different terms in (58). Invoking the character- 
ization of the cones in Lemma 5 and (49), it follows, for example, that for 
(37,22) £ T^2i ^2(x2,T) ~ C2e~ l2X2 whereas for (21,22) 6 T> x D T>i it holds 
that 

(61) ^\x 2 ,T) ~ \D^\v a )\(vJx 2 ) 1 / 2 e-*>#~^-V v ' i 

where D 2 71 ^ is specified by D in Theorem 1 with k = k 2 71 ^ ■ Thus, by 
Lemma 6(i), it follows that the leading term is C2e~" /2X2 if (21,22) G T>2- 
Similarly, it can be checked that the leading term is C\e~^ lXl if (21, 22) € T>\. 
Finally, in the case that (21,22) €T>q, we note that both terms in (58) are 
of the same order [cf. Lemma 6(H)] . More precisely, 

^(zi^M^W + c^^ 

where D 2 , D 2 71 are specified by D in Theorem 1 with k = K2 and k = 

k 2 71 , respectively, and C\ is given in Proposition 2. Using Lemma 6 and 
that k'{ = k 2 , it is a matter of algebra to verify the form of the constants D 2 
and Df. 

Drawing on Proposition 3, Theorem 6 can be proved following an analo- 
gous line of reasoning. We omit the details. □ 

6. Examples. We now develop two explicit examples that illustrate the 
results shown in the previous sections. 

6.1. Cramer-Lundberg model with exponential jumps. Let X be a drift p 
minus a compound Poisson process with rate A and exponential jump sizes 
with mean [i starting at 2. Then the characteristic function of X reads as 
k(9) = p9 — X9 / (fi + 9) and, if p > ^, the ultimate ruin probability is equal 
to tp(x) = Ce~ lx , where the adjustment coefficient is 7 = \x — X/p and C = 
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A/ (up). More generally, it was shown by Asmussen (1984), Knessl and Peters 
(1994) (with p = 1) and Pervozvansky (1998) that the finite time ruin prob- 
ability ip(x,t) is given by 

(62) tf(x,t) = 1 - il>(x,t) = [1 - Ce-^]1 (7>0) + w(x,t), 
where 

(63) w(x,t) = -J— f S+ e a ^ x ~ qt sm(b{q)x - (f>(q)) dq 



vr \j fip J s _ q 
with s± = {y/X±y/JIp) 2 , 4>(q) = arccos(^±=2) and 



(64) a(q)= X -W-\ b(g) = VMtE^E^EW 

Further, we note that, under p(°\ X is still a drift p minus a compound 
Poisson process with exponential jumps with the changed rates A c = A^; 
and p c = p + c. In particular, A_ 7 = pp and p_~ f = X/p are the parameters 
under p( _7 \ 

In view of the previous paragraph, we see that, under P^tO, the drift 
of X2 is always negative, k 2 71 (0) = k' 2 (— 71) < 0. Also, under p(~^ 2 \ the 
adjustment parameter of A"i is positive if and only if p > p* := p\jp\ and is 
then equal to 

(rv\ ~ V ( V\ 

(65) 7 = 73 - 72 = — \p 

P2 \ Pi 

and the asymptotic constant C satisfies C2 = CC2 = Inserting the ex- 
pressions (62)-(64) (with the proper choices of parameters) into Corollary 1 
leads then to explicit expressions for VWL^sim and tp OT . 
It is a matter of calculus to verify that 

pl/p-pi _{pi/p-pi)+ 



Sl = -n-, , S 2 



v\Ip-P2 p%/p-P2 

and, if p > p* , 

_ pp\Ip\-p\ 

*3 2/2 

PPl/P2 -P2 

Also, by invoking Corollary 1 or by a direct calculation, we see that 

Ci(v) = — = Ci, z = 1,2. 

m 

Inserting these quantities into Propositions 2 and 3 yields explicit asymp- 
totics expansions for i/wii V'sim and tp or . 
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6.2. Brownian motion with drift. If X(t) = mt + B(t) where B(t) is stan- 
dard Brownian motion, then its characteristic exponent reads as k(9) = 
^9 2 + mO. If m > 0, ip(x) = e~ yx , where 7 = 2m is the adjustment coeffi- 
cient. Further, under P^ c \ X is still Brownian motion, but the drift changes 
to m + c. The drift of the measure associated to c = —7 is — m, that is, 
the Brownian motion switches its drift. In view of the Corollary 1 and the 
well-known first-passage distribution of Brownian motion with drift, 

X + mt\ _2mx^(~ x + m t 



(66) ^( M ) = $(_^j_ e -*™<^ 

we find that if x 2 > x\ then 

^orOri,x 2 ) = P(n( Xl ) <T)+ e - 2 ^ 2 p(- 2 ^)(r 1 (x 1 ) > T) 
= 1 - $(o(xi,pi)) + e" 2pia;i $(a(-x 1 ,pi)) 

(67) 

+ e~ 2p2X2 x [*(a(x 1 ,p 1 -2p 2 )) 

- e~ 2:El(pi - 2p2) $(a(-xi,pi - 2p 2 ))], 
^ sim (xi , x 2 ) = ^(r 2 (x 2 ) < T) + e" 2 f ™ P^ > (r 2 (xi) > T) 
= 1 - $(a(x 2 ,P2)) + e" 2p2a;2 $(o(-x 2 ,p 2 )) 

(68) 

+ e -2pi*i x [*(o(s 2 ,P2-2pi)) 

_ e -2x 2 ( P2 -2 Pl )$ (a( _ X2i:p2 _2p 1 ))], 

V'and^i,x 2 ) = $(a(x 1 ,p 1 ))-l + e- 2pia;i (l-^(«(-^i,Pi))) 

(69) + e~ 2p2X2 x [e" 2a:i(pi - 2p2) $(a(-x 1 ,p 1 - 2p 2 )) 

+ l-$(a(x 1 ,pi-2p 2 ))], 

where a(x,p) = [x+j>T]/\/T and $ denotes the cumulative standard normal 
distribution function. In view of the facts that 3>(— x) = 1 — $(x) and 1 — 
$(x) ~ (2-7t)~ 1 / 2 x _1 exp(— x 2 /2) as x — > 00, it follows from (68) and (69) that 
if xi,x 2 tend to infinity with x\jx 2 = a then 

e -2 P2 x 2 -2(p 1 -2 P2 )+ a;i + ( 0Va ) ) if < a < s 3 , 

(70) ipandixi, x 2 ) « ^ o(av a ), ifs 3 <a<s 1 , 

e -2pixi _|_ ( aUa ) ) ifsi<a<l, 

e -2w*a+(y(t; a ), if0<a<s 2 , 

(71) Aim(xi,x 2 ) ^ { 5(v a ), ifs 2 <a<si, 

e -2pi*i +J( Uo ), if Sl <a<l, 

where 
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S(v) 



2v 



+ 



2v 



y2 ~ P2 (P2 ~~ 2pi) 2 — V 2 



and 



si 



Pi 



2pi - p 2 ' 



S2 



\/27rx2 
{2p 2 -pi) 

P2 



-x 2 {v+p 2 ) 2 /[2v] 



and, if pi > 2p2, 



S3 



Pi - 2p 2 



2pi - 3p 2 

The asymptotics of V'sim agree with the asymptotics of the steady state 
distribution of a tandem queue calculated in Lieshout and Mandjes (2007). 

By straightforward calculations, it can be verified that if S is a Brownian 
motion, then Cj = 1 and 8 V = —v — P2, 0' v = v — P2, 9* = v + p2 — 2p\ and 



Ciiy) 



1. 



«?(-«) 



>n(-at>a) = /^(-Va) 



for 2 = 1,2. Inserting these quantities in (58)— (59) and comparing with (70) 
and (71) it follows that Propositions 2 and 3 and Theorems 5 and 6 remain 
valid if S is a Brownian motion. 



APPENDIX 

In this section, we briefly review the general large deviations theory for 
first passage times developed by Collamore (1996) and explicitly relate it to 
the results for i/'sim in Theorem 5 by calculating the relevant quantities of 
the general theory. 

We denote by (X\, X2) the two-dimensional Levy process given by Xi(t) = 
Pit — S(t) where S is a Levy process (starting at zero) and write A = 
A a = (— 00,— a) x (— 00, — 1) and Tg for the first hitting time of a set B 
by (X\, X2). The probability tp s i m (aK, K) is then equal to P{Tka < 00). If 
^($1)02) = logE[e dlXl+8 ' 2X2 ] is finite in a neighborhood of the Cramer set 
C = {6 : k(0) < 0} and G C°, then it holds that [Theorem 2.1 of Collamore 
(1996)] 

Jim —\ogtjj B i ra (aK,K) = - inf I(x), 

K-*oo K xeA a 

where / is the support function of the Cramer set C, 
(72) J(x)=sup(0,x). 

In the next result I(x) is identified: 
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Proposition 4. For x\,x 2 < 0, it holds that 

I(x 1 ,x 2 ) = \x 2 \j(xi/x 2 ), 

where 



7(a) 



■inf(0,l) 



ifa^l, 
ifa=l 



with 1 = (1,1), k 2 (s) := supggjjfls - n 2 (0)} and v a = (p 1 -p 2 )/(l - a). 
Using this result, we can calculate the first-passage rate function: 
Corollary 4. It holds for a > that 



72, if0<a<s 2 , 
inf I(x) = { 7(a), if s 2 < a < s\, 
«7i, ifa>si. 



X&A a 



We observe here that Corollary 4 agrees with the exponents of the asymp- 
totics found in Theorem 5 (as it should). 

Proof of Proposition 4. The linear functional 9 i— ► (9,x) attains it 
maximum over the closed set C at a point 9* of the boundary dC. In view of 
the form of k (Example 1), the definition of v a and k(9*) = 0, it is a matter 
of algebra to check that for a / 1 

V„ 



I(xi,x 2 ) = sup< Q\xi + 6» 2 x 2 - («2(0i + #2) + (pi -pz)0i] 
e&ac 1 



£2 

»2 
On 



sup{(ei + e 2 )(-KI)-«2(e 1 + 02)} 

<9e<9C 



sup{7?(-|u |) - K 2 (r])} 
r{>e_ 



;t2_ 



M), 



where 9 = ini{9 : k 2 (9) < 00} and in the last line we used that k 2 (9i + 9 2 ) = 
—01 (pi — P2) for (9\,9 2 ) G dC. The rest of the statements follows by straight- 
forward calculations. □ 



Proof of Corollary 4. From the form of T, we deduce that / attains 
it minimum over A a at the boundary dA a . If we set x 2 = —1, the minimiza- 
tion reduces to inffe> a 7(6). Taking note of the fact that k 2 (—v) > 72^ with 
equality if and only if v = — s 2 we see that 



(73) 



inf 7(6) 

b>a 



72, 

7(a), 



if < a < s 2 , 
if a > s 2 . 
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Similarly, setting x\ = —a, leads to the minimization inf c <i cj(a/c) or equiv- 
alently, inf<2<<j7(c2)/d. Observing that 



we see that k^ - \ v a\) > 07iba| with equality if and only if a = si and con- 
clude that 



Combining equations (73) and (74) completes the proofs. □ 

Let us also note that classical arguments for the large deviations (LD) 
theory for stationary increments processes allow us to explain heuristically 
the structure of "or" and "sim" ruins in the quadrant. Let us consider for 
example the "sim" ruin. From the LD theory, we expect that paths exiting 
the positive quadrant, seen from far away, will be concentrated near one of 
three possible directions: the direct path to the origin, and the "dominant" 
(most probable) paths reaching the x\ = 0/x2 = axes, respectively. 

It turns out that all the dominant exit paths for reaching the {x± < 0, x<i = 
0} semi-axis are parallel to each other and, therefore, give rise to a "bound- 
ary" cone T>\. Similarly, the dominant exit paths reaching the x\ = 0,X2 < 
semi-axis gives rise to a cone T>2- These two cones will be disjoint in our 
case, as indicated in Figure 2. In these two boundary cones, the probability 
of simultaneous ruin is equivalent to the probability of ruin of the X1/X2 
process respectively, that is, it holds asymptotically for large x±, X2 on a ray 
that lies within these cones T>i that V'sim^i, £2) ~ ipi(%i)- 

A similar result holds for the "or" ruin. We may note geometrically that 
the boundary cones for hitting the semi-axes {x\ = 0, X2 > 0} and X2 = 
0, x\ > are precisely the complements C2 = T>\, C\ = V c 2 . 

This intuitive picture is confirmed and sharpened by Theorem 5 and 
Proposition 2. 
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